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r^ Abstract: Diffusion coefficient measurements are important for many bi- 

K- 5 ological and material investigations, such as particle dynamics, kinetics, and 

T-H size determinations. Amongst current measurement methods, single particle 

^^ tracking (SPT) offers the unique capability of providing location and diffu- 

I— I sion information of a molecule simultaneously while using only femptomoles 

^ of sample. However, the temporal resolution of SPT is limited to seconds 

PQ for single-color labeled samples. By directly imaging three dimensional (3D) 

Q diffusing fluorescent proteins and studying the widths of their intensity pro- 

' r^ files, we determine the proteins' diffusion coefficients using single protein 
images of sub-millisecond exposure times. This simple method improves the 
temporal resolution of diffusion coefficient measurements to sub-millisecond, 

,__( and can be readily applied to a range of particle sizes in SPT investigations 

P' and applications where diffusion coefficient measurements are needed, such 

P^ as reaction kinetics and particle size determinations. 

Determination of particles' diffusion coefficients is important for many bi- 
ological and material applications, such as single-molecule dynamics studies 
O [ini El |9j, biolochemical and pharmaceutical reaction kinetics studies [H [5], 

and particle size and shape determinations [6]. Amongst current methods for 
;L| measuring diffusion coefficients, which include nuclear magnetic resonance 

•^ imaging (NMR) [7j, dynamic light scattering (DLS) [B], fluorescence corre- 

rN lation spectroscopy (PCS) PUnilH], and fluorescence recovery after photo- 

a bleaching (FRAP) [12], single particle tracking (SPT) offers the unique capa- 

bility of simultaneous location and diffusion coefficient determination. This 
is essential for molecular mechanism investigations in heterogeneous environ- 
ments, such as inside a cell's cytoplasm [13j, flagella [H], and membrane [U] 
in vivo, and on DNA molecules [lO] in vitro. Due to this capability, and the 
additional advantage that SPT experiments require less than femtomoles of 
sample, SPT can be a powerful tool in diffusion coefficient measurements 
applicable to a large range of biological investigations in vitro and in vivo 



where supplies are scarce. 

However, the drawback of using SPT for diffusion coefficient measure- 
ments is the low temporal resolution. In single-molecule fluorescence imag- 
ing studies, stationary or slowly moving (relative to the data-acquisition 
timescales) single molecule intensity profiles are called point spread func- 
tions (PSF), and are fit to Gaussian functions to determine the molecules' 
localization information. The centroid of the Gaussian function determines 
the lateral location of the molecule at the time of imaging, and the standard 
deviation (SD) determines the axial location. In SPT diffusion coefficient 
measurements, consecutive locations of a single fiuorophore are measured and 
from mean square displacement analysis of the particle's single trajectories, 
diffusion coefficients are obtained [101 [13], [16] . This method requires at least 
20 consecutive location measurements for each single trajectory. With the 
current single-photon camera imaging rate of approximately 100 frames/sec 
for a finite-sized imaging area, 0.2 sec is required and 3D diffusion coeffi- 
cient D^ measurements up to order 10^ nm^/s have been reported [T7|. This 
requirement of 0.2 sec is, however, too long for diffusion coefficient measure- 
ments of fast- moving molecules, such as nanometer-sized proteins that diffuse 
beyond the typical imaging depth of ~ 400 nm of single-molecule imaging 
microscope setups in less than 1 ms (a typical 5 nm protein has D^ ^ 10^ 
nm^/s and diffuses yj2D^t ~ 447 nm in 1 ms). A recent SPT method mea- 
sures D^ up to 1.7 xlO^ nm^/s by labeling the particles with two colors 
[18]; however, multi-color-labeling may not be feasible for many biological 
particles of interest and thus restricts the applicability of the method. 

A SPT method that can determine 3D diffusion coefficients of single- 
colored nanometer-sized biological entities in their native environment is 
highly desirable for in vivo and in vitro studies. In order to capture the 
molecule within the microscope's imaging depth, the imaging time needs to 
be less than 1 ms. Here we report a novel method that determines the diffu- 
sion coefficient of nanometer-sized Brownian molecules from the SD values of 
the molecules' intensity profiles using sub-millisecond exposure times. This 
is a single-image molecular analysis (SIMA) study of dynamic molecules and 
is an extension of our previous stationary molecule investigations [3]. In 
this study we used eGFP as the nanometer-sized fluorescent molecule for 
measurements and analyses. 

Since the imaging times in our method are less than 1 ms, the temporal 
resolution of diffusion coefficient measurements is improved by at least 1000- 
fold over the minutes-long PCS method (multiple measurements each of order 



20-second long), 200-fold over the 0.2-sec long centroid SPT method, 50- 
fold over the typically 50-nis long FRAP method, and 10-fold over the two- 
color SPT method. Furthermore, the improvement in temporal resolution 
is achieved without compromising the precision of the D^ measurements, 
and the single-image nature of the method avoids the photobleaching and 
limited lifetime photon problems associated with single-molecule fluorescence 
imaging studies. Below we show our measurement method that relates the SD 
of a 3D freely diffusing protein's intensity profile to its diffusion coefficient 
1^3. A prior study has used a similar concept to relate slow 2D diffusion 
coefficients (up to 1.1 x 10^ nm^/s) to a fiuorophore's spot sizes [20]; here 
we extend the study to fast 3D diffusion. 

Figures [l] and [2] illustrate the principle of this method. In a finite ex- 
posure time, the intensity profile of a moving molecule is wider (or more 
blurry) compared to that of an immobile molecule. Figure [IK shows a 30 ms 
frame image of stationary eGFP molecules adsorbed on a fused-silica surface, 
and Fig. [Tp shows a 1 ms frame image of diffusing eGFP molecules near a 
hydrophilic fused-silica surface [1]. These figures clearly show that diffusing 
molecules images are blurry compared to that of immobile molecules. In 
Figs. [lj3 and [lj5 the intensity profiles of the stationary and diffusing eGFP 
molecules are plotted, and in Figs. [Tp and [Tp, the respectively selected in- 
tensity profiles are fitted to Gaussian functions. While both intensity profiles 
fit well to a Gaussian function, the width (or SD) of the diffusing protein's 
intensity profile is larger than that of the stationary protein. 

In general, the final image of a diffusing molecule, like those in Fig. [1} 
is the sum of the emitted photons along its diffusion trajectory projected 
onto a 2D imaging screen during the exposure time. Figure [2|\ shows a 
simulated eGFP diffusion trajectory at 0.6 ms exposure time using 0.005 ms 
steps for clarity. The data are gray-scaled to correspond to the particle's 
axial locations (SOM text). The emitted photons, after photon-to-camera 
count conversion, were projected onto a 2D imaging screen and binned into 
our camera pixels each of 79 x 79 nm^ in size (Fig. pp; gray image) and 
the corresponding diffusing eGFP PSF intensity profile was formed in the 
colored image above. The total photon count of this image was 414. The 
2D Gaussian fit to the diffusing eGFP intensity profile is shown in Fig. |2p, 
yielding SD values in the x- and ^/-directions, Sx and Sy, respectively. Sx,y 
values presented in this article are results from fitting to these experimental 
and simulated PSF data, and they were used to quantify the blur of diffusing 
eGFP molecules, and consequently, the diffusion coefficient D^. 



In order to determine D3 from diffusing fluorophore images, we performed 
(i) experimental measurements, (ii) analytical calculations, and (iii) simula- 
tions and obtained SDs of diffusing eGFP intensity profiles at 0.3 to 1 ms 
exposure times. Below we show that when the experimental results were 
checked against theoretical calculation and numerical simulation results in 
Fig. |4]A., the good agreement validates our method of measuring nanometer- 
sized fluorophore diffusion coefficients. 

In experimental measurements. Fig. [3]\ shows representative eGFP im- 
ages (chosen such that the molecule's respective Sx values were within ± 5 
nm of the means to the respective diffusing eGFP intensity profile SD distri- 
butions in Fig. IsIB) at 0.3, 0.7, and 1 ms exposure times. As expected, SD 
values of these respective single diffusing eGFP molecules increase from 136.4 
to 160.9 and 175.5 nm, validating that SD provides a quantitive measure for 
the motion-induced blurriness of single fluorophore images. The increasing 
mean eGFP SD values with exposure time are plotted in Fig. |4]A., where the 
error bars are the SDs of the eGFP intensity profile SD distributions in Fig. 

m. 

In analytical calculations, we deduce an expression relating a diffusing 
eGFP's SD to D^. The study involves first projecting the eGFP PSFs at all 
focal depth onto a 2D imaging screen, forming an axial-direction-projected 
PSF /(x, y), and then convolving this projected PSF with the lateral location 
distribution of the molecule in a trajectory, which we define as a pathway 
distribution function (PWDFj^y for the lateral directions) g{x,y): 

I{.x,y)(x.f{x,y)*g{x,y). (1) 

In SOM we show that both /(x, y) and g{x, y) can be approximated well 
by Gaussian function for sufficiently low exposure times, and their convolu- 
tion is another Gaussian function with a variance equal to the sum of the two 
variances. Therefore, the final projected intensity profile of a 3D diffusing 
molecule is a Gaussian function with SD being 



Sx^y = \ s'^ + Ax,y ■ 2D3t, (2) 



where s'o(t) = ^/lll^TOMMD^ ^ ^/sl + OMSiDst nm is the SD of the 
axial-direction- projected PSFs for our experimental parameters, and A^^y ■ 
2D^t is the variance of PWDFj, j^s with A^^y = 0.0926. This relation enables 
the determination of D3 from the SD of a single-molecule's intensity profile 



and the exposure time as 



s 2 



D = ^'^ (3) 

^ (2v4^,y + 0.0634)t' ^ ^ 



Using FCS-determined eGFP D3 = 8.86 x 10'' nm^/s in Eq. [8| the ana- 
lytical eGFP Sx results are plotted in Fig. |4]A., showing excellent agreement 
with the experimental Sx results within 0.7 ms. Note that Sx starts to devi- 
ate from the experimental results at t > 0.8 ms; this is because the exposure 
time begins to approach the diffraction-limit-determined value for eGFP for 
this study (SOM text). 

We have also performed simulations of diffusing eGFP intensity profiles 
(as shown in Fig. [2]) using the FCS-determined D^. Figure EK juxtaposes the 
simulated diffusing eGFP SD results with the experimental results; the two 
mean values and error bars agree at all exposure times (Fig. S7 compares 
the results at t = 0.6 ms). 

To determine the precision of the measured D^ from single eGFP images, 
we performed error propagation analysis of D^{sx,y) using Eq. [s] 
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ADg = — ^ ^ASxv, (4) 

^ {Ax,y + 0.032)t '^'^' ^ ^ 

where Asx,y is the SD measurement precision of the single fluorophore's in- 
tensity profile (the experimental error bars in Fig. IlK, also Ref. [2]). Figure 
|4p compares the experimentally determined D^ and AD^ from single diffus- 
ing eGFP image SD measurements to the FCS-determined eGFP D^ = 8.86 
nm^/s, showing agreement. 

At 0.7 ms, AD3 = 5.2 xlO^ nm^/s for a single eGFP image using both 
the statistically independent mean Sx and Sy values of {sx,y) = 162.1 nm 
and Asx,y = 39.2 nm. It is 57% of the eGFP D3 value of 8.86 xlO^ nmVs. 
Since there are order 30 molecules in a typical frame image of less than 1 
ms exposure time, the precision to D^ measurement further improves by 
\/30 times to 10%, comparable to the precision of FCS D3 measurements 
[5]. In spatially restrictive situations, such as in vivo imaging in typically 
micron-sized cells, where only one image can be obtained at a time, repeated 
single-image measurements will enable precise D3 determination. 

Although this study focuses on fast diffusion of nanometer-sized proteins 
in free solution with D3 > 5 x 10^ nm^/s, the methodology applies to 3D dif- 
fusion of all rates. When diffusion coefficients are low for large particles, in a 



crowded environment, or in viscous solvents (such as in cells [23] or glycerol), 
the molecule's intensity profile will be more localized. Consequently, longer 
exposure times should be used to observe noticeable changes to the SD from 
the stationary values. SOM text explains the procedure in determining the 
appropriate exposure times for a particle of unknown D^. 

Molecules' movement can deviate from a 3D unbiased Brownian motion, 
examples include directional motion and diffusion with drift. For these al- 
ternative motions, s'q and PWDF^^^^ ^ should be determined before convolv- 
ing the axial-direction-projected PSF with PWDFa;^^^ for the final intensity 
profile. As long as the mean numerical SD of locations in the molecule's 
trajectory is less than half of the diffraction limit at the exposure time, the 
projected convolved image of the molecule will be a unimodal intensity profile 
that can be fitted to a Gaussian function, and the resulting Sx,y will provide 
information on the molecule's dynamics. 

In summary, we present a new single-molecule fluorescence image analy- 
sis method that measures fast diffusion coefficients with high precision. The 
experimental setup and data analysis are simple while using standard mi- 
croscopy imaging systems, and the method is applicable to a wide range of 
diffusion coefficient measurements with greatly improved temporal resolu- 
tion. Applications in basic research and pharmaceutical investigations such 
as fast drug screening can be envisioned. 
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Figure 1: Comparing stationary to diffusing eGFP molecules. (A) An im- 
age of stationary eGFP molecules adsorbed on a fused-silica surface. Five 
of the seven molecules have signal-to-noise ratio (SNR) > 2.5. (B) Intensity 
profiles of the stationary eGFP molecules in (A) in photon counts. (C) Inten- 
sity profile (dots) and Gaussian fit (mesh) to the stationary eGFP molecule 
denoted by arrow in (A) and (B). For this molecule, the SNR is 9.8, Sx = 
107.2 nm, and Sy = 107.9 nm. (D) Diffusing eGFP molecules near a reflec- 
tive hydrophilic fused-silica surface at 1 ms exposure time. Six of the eight 
molecules have a SNR > 2.5. The scale bars for (A) and (D) are 2 /im. (E) 
Intensity profiles of the diffusing eGFP molecules in (D). (F) Intensity profile 
(dots) and Gaussian fit (mesh) to the diffusing eGFP molecule denoted by 
arrow in (D) and (E). For this molecule, the SNR is 3.5, Sx = 202.2 nm, and 
Sy = 192.4 nm. It is clear that the intensity profiles of diffusing molecules 
are wider (or have larger SDs) than that of stationary molecules. 
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Figure 2: Simulated image formation and analysis process of a diffusing 
eGFP molecule. (A) Trajectory of a diffusing eGFP molecule in free solution 
under TIRF (total internal reflection fluorescence) evanescent excitation at 
the exposure time of 0.6 ms. The data is gray-scaled to correspond to the 
particle's axial locations (SOM Text). (B) The emitted photons from the 
trajectory form an intensity profile (colored plot), which is then projected 
onto a 2D camera screen (black and white image). (C) Gaussian fit (mesh) 
to the intensity profile of the diffusing eGFP (dots), where s^ = 119.4 nm, 
and s„ = 142.2 nm. 
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Figure 3: Diffusing eGFP images and intensity profile SD distributions at 
different exposure times. (A) Three representative images showing diffusing 
eGFP molecules at exposure times of 0.3, 0.7, and 1 ms. The intensity profile 
SD values increase with the exposure time. The scale bar is 1 yum. (B) EGFP 
intensity profile SD distributions (normalized by counts for comparison) at 
the three aforementioned exposure times, showing increasing values of 136.8 
± 27.7 (mean ± SD), 159.0 ± 32.24, and 172.1 ± 34.8 nm, respectively. 
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Figure 4: Comparing Sx and D^ results. (A) Experimental (circles), simula- 
tion (disks), and theoretical calculation (squares) measurements of diffusing 
eGFP intensity profiles' mean Sx vs t. In the experimental and simulation 
results, the error bars are the SDs of the Sx distributions. (B) Experimental 
D3 calculated from Eq. [3] The error bars are AD^ calculated using Eq. 
|4J the dashed hue is the FCS-determined eGFP D^ of 8.86 xlO'' nm^/s for 
comparison. 
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Supporting Online Material 
Materials and methods 

Sample preparation and imaging 

The eGFP molecules (4999-100, Bio Vision, Mountain View, CA) were di- 
luted in 0.5X TBE buffer (45 mM Tris, 45 mM Boric Acid, 1 mM EDTA, pH 
8.0) to 0.03 nM. For the stationary eGFP studies, manufacturer pre-cleaned 
fused-silica chips (6W675-575 20C, Hoya Corporation USA, San Jose, CA) 
were used, where isolated eGFP molecules were adsorbed to surfaces at low 
concentration. For the diffusing eGFP studies, the manufacturer pre-cleaned 
fused-silica chips were treated using oxygen plasma for three minutes, ren- 
dering it hydrophilic to prevent eGFP adsorption [T]. The hydrophilic fused- 
silica surface can be considered ballistic for the diffusing eGFP molecules in 
our experiments and simulations. For both studies, a protein solution of 5 
/iL was sandwiched between the fused-silica surface and an oxygen-plasma- 
cleaned coverslip (2.2 x 2.2 cm^), resulting in a 10.5 /xm thick water layer. 
The coverslip edges were then sealed with nail polish. 

Single-molecule imaging was performed using a Nikon Eclipse TE2000- 
S inverted microscope (Nikon, Melville, NY) in combination with a Nikon 
lOOX objective (Nikon, 1.49 N.A., oil immersion). Samples were excited by 
prism-type total internal reflection fluorescence (TIRF) microscopy with a 
linearly polarized 488 nm laser line (I70C-SPECTRUM Argon/Krypton laser. 
Coherent Inc., Santa Clara, CA) focused on a 40 x 20 /im^ region. The 488 
nm line was filtered from the multiline laser emission using a polychromatic 
acousto-optic filters (48062 PCAOM model, NEOS Technologies, Melbourne, 
FL). The laser excitation was pulsed with illumination interval of 30 ms for 
the stationary eGFP molecules in Fig. 1 and Fig. ^ and between 0.3 
ms and 1 ms for the diffusing eGFP molecules. The excitation intensities 
were 2.7 and 3.2 kW/cm^ for the respective stationary eGFP molecules, 
and 37.5 kW/cm^ for the diffusing molecules. Images were captured by an 
iXon back-illuminated electron multiplying charge coupled device (EMCCD) 
camera (DV897ECS-BV, Andor Technology, Belfast, Northern Ireland). An 
additional 2X expansion lens was placed before the EMCCD, producing a 
pixel size of 79 nm. The excitation filter was 488 nm/10 nm, and the emission 
filter was 525 nm/50 nm. 

Data acquisition and selection 

Movies were obtained by synchronizing the onset of camera exposure with 

laser illumination for different intervals. The maximum gain level of the 
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camera was used and the data acquisition rate was 1 MHz pixels/sec (~ 3.3 
frames/sec). Single-molecule images were checked such that there were no 
saturations in the intensity profiles. For the defocusing analysis of stationary 
eGFP molecules, 21 x 21 pixel boxes centered at the molecule were selected 
by hand using Image J (NIH, Bethesda, MD), and the intensity values were 
used for the 2D Gaussian fitting. For the diffusing eGFP molecule movies, all 
visible diffusing eGFP intensity profiles in the peak laser excitation region of 
10 X 10 /im^ were selected by hand using 39 x 39 pixel boxes centered at the 
molecule. The center 25 x 25 pixels of the boxes were used for 2D Gaussian 
fitting, and the peripheral pixels were used for experimental background 
analysis. 

Before analysis, the camera's intensity count at each pixel in an image 
was converted into photon count by using the camera-to-photon count con- 
version factor calibrated the same day of the measurement as described in 
our previous article p]. The number of detected photons in an image was 
obtained by subtracting the total photon count of the EG from the total 
photon count of the image. The eGFP intensity profiles were fit to a 2D 
Gaussian function in order to obtain the SD values of the molecule: 



f{x,y) = /oexp 



X - xof (y - VoY 



2sl 2sl 



+ (&), (5) 



where /o is the multiphcation factor, s^ and Sy are SDs in the x- and y- 
directions, respectively, xo and yo are the centroid location of the molecule, 
and (b) is the mean background offset in photons. 

For the defocusing eGFP analysis, we selected 17 adsorbed eGFP molecules 



with a minimum photon count of 229 and signal-to-noise ratios (SNR, Jo/y /q + cr^) 
of higher than 3.75, where Iq is the peak PSF photon count (after subtract- 
ing the mean background offset (b)) and a^ is the background variance in 
photons. For the diffusing eGFP molecules, we used a SNR of 2.5 as a se- 
lection criterion, and PSFs with photon counts less than 50 were not used 
in the analysis because they would be invisible in experimental data. At 
each exposure time, we acquired 1600 data points from 4 movies recorded 
at different regions of the imaging chip; the number of diffusing eGFP data 
used for the experimental analysis that satisfied the SNR criteria are 419 to 
1066 for the 0.3 ms to 1 ms exposure times, respectively. 



^x) 



Analytical expression of diffusing eGFP D3,{sa 

In this section we decompose an eGFP's 3D diffusion process into two com- 
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ponents for s^ and D^ calculation: a ID diffusion along the axial direction 
and a 2D diffusion in the lateral directions. 

It is known that as the defocusing distance between the fluorophore and 
the focal plane increases the PSF's SD increases as well. Consequently, cal- 
culation of the intensity profile necessitates integrating over all axial loca- 
tions the molecule may have traveled during the exposure time to obtain an 
axial-direction-projected PSF, f{x,y). As diffusion in the lateral and axial 
dimensions are statistically independent of each other, we choose to perform 
this integration prior to convolving the resulting PSF with PWDF^^^y in the 
lateral dimensions to obtain the final projected 2D intensity profile of the 3D 
diffusing molecule on an imaging screen. 

In the axial direction, the axial-direction-projected PSF is computed by 
numerically integrating defocused PSFs through z for all pixelated x, y values 
(sufficient to contain all defocused PSFs) 



400 

C{z) exp ( — 








dz. 



2s,{z)^ 2sy{zy 

(6) 
where C{z) and Sx,y{z) are the amplitude and SDs of our imaged, defo- 
cused eGFP Gaussian PSFs (SOM text), respectively, (2:0) and A^ ■ 2D^t 
are the mean and variance of diffusing eCFPs' Gaussian PWDF^s (SOM 
text), exp{—z/d) describes the decaying TIRF evanescent excitation inten- 
sity, and the range for the 2;-integration is the imaging depth of nm to 400 
nm measured from the focal point at the fused-silica surface. The result- 
ing axial-direction-projected PSF f{x,y) remains Gaussian, and the SD is a 
function of the exposure time t as SQ(t) = i/lll^ + 0.0634_D3t nm. 

In the lateral directions, the intensity profile I{x, y) of a diffusing eGFP's 
image is the convolution of the axial-direction-projected eGFP PSF f{x, y) 
with the PWDFa;^j^, g{x,y) as 

I{.x,y)(x.f{x,y)*g{x,y). (7) 

We numerically calculate g{x,y) of a freely diffusing eGFP particle by 
simulations. Figure S|5]A. shows 9 random PWDF^jS at exposure time t = 0.6 
ms. Six of the nine PWDF^jS have one peak (uni-peaked or unimodal) and 
can be fitted to a Gaussian function with R'^ > 0.8. Figure qSp shows the 
SD distribution of PWDF^-s, combining the Gaussian fitted SD values for 
the uni-peaked PWDF^s and the numerical particle location distribution SD 
values for the double-peaked PWDF^jS; the mean is 96.8 nm. Figure ^5p 
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shows that when the 9 PWDF^-s in Fig. ^5]A. are convolved with single- 
eGFP PSFs at focus with sq = 108.2 nm, all convolved PWDF^s fit well 
to a Gaussian function, and the mean of the SD distribution is 147.1 nm. 
Therefore, although not all PWDF^.s are uni-peaked, taken over all, we can 
view PWDFa-s as Gaussian functions with an average t-dependent SD value 
of y/Ax-2D:i,t. For the 0.6 ms exposure time data, A^ = 0.0882. We found 
Ax to be insensitive to exposure times below 1 ms (mean A^ is 0.0926; data 
not shown). 

Given that /(x, y) (at focus and the axial-direction-projected) and g{x, y)s 
are both Gaussian functions, in the lateral directions, their convolution can 
be described by another Gaussian function with a variance equal to the sum 
of the two variances. Using the focused eGFP PSFs with sq = 108.2 nm and 
PWDF^ at 0.6 ms, 3^,20 = ^/sl + A^ ■ 2D^t = ^108.22 + 96.8^ nm = 145.2 
nm, very close to the mean SD value of the above PSF-convolved-PWDF^s 
of 147.1 nm. 

Finally, we can calculate the intensity profiles' SDs of freely 3D diffusing 
molecules. Since we have observed that both the axial-direction-projected 
PSFs (data not shown) and the lateral PWDFi.^j^s are Gaussian, the final 
projected intensity profiles' SD of diffusing molecules is 



Scc,y = V ^0' + ^-.S/ ■ 2^3^. (8) 

This relation enables the determination of D3 from the SD of a single- 
molecule's intensity profile and the exposure time. 

Diffusing eGFP Simulations 

We simulated 3D Brownian diffusion eGFP trajectories at a range of exposure 
times using FCS-determined eGFP D3 = 8.86 x 10^ nm^/s and triplet state 
statistics. The starting locations of the trajectories followed the distribution 
function described in the SOM Text. The step sizes in the x,y, and z direc- 
tions were randomly selected from a Gaussian distribution with a mean of 
zero and SD of y/2D3to with a step time to = 1 A^s. Because of the reflective 
fused-silica-water interface, the simulated z values were maintained above 
zero. The number of steps in a simulation was t/to- At each x, y location 
in a trajectory, when the molecule was not in a triplet dark state, a Poisson 
distributed number of photons (described in SOM text) were drawn from a 
Gaussian PSF spatial distribution with a mean of zero and the corresponding 
SD value for the axial-location (SOM text). This relative displacement of the 
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photons is added to the simulated x, y location of the molecule, generating 
the actual x, y location of the emitted photons at the simulation step. 

The simulated photons of each trajectory were binned into 50 x 50 pixels 
with a pixel size of 79 nm. Then the photon count of each pixel was converted 
into the modified camera count using Eq. 4 in Ref. [2] with the photon mul- 
tiplication factor of the camera M = 1 in order to include the camera count 
variance effect. Random background photons at each pixel were generated 
using the corresponding experimental background distribution functions for 
the exposure time (described in Ref. [2]). The final intensity profiles were 
fit to a 2D Gaussian function to obtain the two SD values for the image. 
For each SD datum of diffusing eGFP molecules in Fig. 4, 1000 independent 
trajectories were simulated. 



SUPPORTING TEXT 

Exposure time limit determination for D3 measurements using Eq. 
3. 

In Fig. 4, the calculated SD starts to deviate from the experimental and 
simulation results at 0.8 ms. This suggests the existence of an upper bound 
exposure time for our eGFP studies. In order to determine the appropriate 
exposure time, we explain the origin of this deviation at long exposure times. 

The PWDFajS in Fig. S|5]A. show both unimodal and double-peaked pat- 
terns. When convolved with PSFs, at short exposure times, both the uni- 
modal and double-peaked PWDF^-s will produce unimodal intensity profiles 
suitable for Gaussian fitting; however, at long exposure times, the two-peaked 
PWDFa;S will yield a two-peaked intensity profile. As the exposure time 
increases, both the fraction of two-peaked PWDF^jS and the peak separa- 
tions increase, while Eq. [6] still assumes Gaussian fits for all PSF-convolved- 
PWDFa^.s. As a consequence, a deviation between the analytical and the 
experimental s^ appears and increases with exposure time. 

The threshold exposure time for the onset of the deviation is determined 
as follows. When two identical fluorophores are separated by more than the 
diffraction limit, the combined intensity profile is double-peaked ^. How- 
ever, when the concentration of fluorophores peaks in between the diffraction 
limit separation of the instrument, as for unimodal PWDF^-s, the combined 
intensity profile appear to be unimodal, and can be fitted to a Gaussian 
function. For the double-peaked PWDF^jS, which creates two clusters of flu- 
orophores separated by the distance between the two peaks, this distance, 
which is approximately the numerical 2xSD value of the molecule's location 
distribution, or 2 ■ ^/Ax■2D3t of the PWDF^s, determines the exposure time 
limit for when the convolved intensity profiles become two-peaked. When 
2 • \/Ax-2Dst is less the the diffraction limit of 217 nm, the convolved in- 
tensity profile is unimodal and can be approximated by a Gaussian function. 
The threshold is crossed at the exposure time t ~ 0.8 ms, where the PWDFa;'s 
mean SD value is 113.6 nm. 

From the above analysis, for D^ determination using Eq. 3 in the main 
text, the upper bound exposure time can be determined by requiring 2 • 
y/Ax-2D^t to be less than the diffraction limit separation of the imaging 
setup's emission wavelength and N.A. 

Choosing exposure times for a particle of unknown D3. 

For particles of unknown D^, exposure time can be scanned until the diffusing 
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particle images are noticeabily larger than that of stationary particles, while 
remaining unimodal. In this range of exposure times, we can use Eq. 3 
for eGFP to calculate D3 of this particle for the following reason: in Eq. 
3, since Sg is calculated from integration that depends on D^t (SOM), and 
A^^y ■ 2D^t varies with D^t only, at the appropriate exposure time, although 
the D^ values of the particle and eGFP are different, the D^t values can be 
equivalent. At these exposure times, Eq. 3 for eGFP is restored valid for the 
particle and can be used to determine the unknown D^. 

FCS determination of the eGFP diffusion coefficient. 

In order to independently verify our experimentally determined eGFP PSF 
mean SD results (Eq. 3) ; therefore also the D^ result) by using theoretical 
calculations and simulations, we performed FCS -D3 measurements of eGFP 
(at the Washington University Fluorescence Correlation Spectroscopy and 
Confocal Imaging Facility in the Department of Biochemistry and Molecular 
Biophysics). 

In FCS measurements, fluorescence from freely diffusing eGFP molecules 
at 3 nM concentration in 0.5X TBE buffer (pH 8) was measured. An auto- 
correlation function was used to obtain the eGFP diffusion parameters [1], 

1 [l-F + Fexp{-^)\ 
G{t) = , , p " +1, (9) 

where r is the detection time, A^ is the number of molecules in the detection 
radius w, s is the structure parameter of the excitation beam focal region (the 
ratio of the beam radius in z to the beam radius in x and y), t^ = -fjy is the 
molecule's diffusion time in the imaging area, F is the fraction of molecules 
in triplet state, and r^ is the triplet state lifetime. 

The excitation wavelength for the FCS measurement was 488 nm, and the 
emission photons went through a 505 - 550 nm filter. The excitation power 
was 76.4 kW/cm^, which was comparable to our excitation power of 37.5 
kW/cm^ in the diffusing eGFP studies. We used Alexa 488 with a known 
diffusion coefficient D^ = 4.35 x 10® nm^/s j5j for calibration and obtained 
w ~ 250 nm. For Alexa 488, Td = 35.6 /is. Figure S2 shows the G{t) vs. 
t plot of the eGFP system, where r^ was 174.8 /xs. Assuming a Gaussian 
detection volume and using a one-component fit, the best fit to Fig. S [9] 

2 — 

yields F = 12.7, r^ = 3 /xs, and s = 10. Using r^ = fjy, we obtained eGFP 
D3 = 8.86 X 10^ nm^/s. This value of eGFP D^ is consistent with reported 
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values [5]. 

Refractive index mismatch corrections to eGFP intensity profiles. 

When a fluorescent molecule in water is imaged through a glass coverslip 
using a high N.A. oil immersion objective, the refractive index mismatch 
between the water-based solvent and the glass coverslip changes the fluo- 
rophore's intensity profiles in two major ways [6]: (i) Due to Snell's law of 
refraction, the actual axial location of the molecule (measured from the glass 
coverslip-water interface) is deeper than the apparent axial position of the 
molecule (defined by the depth in water where the imaged fiuorophore's PSF 
amplitude is maximal), (ii) Due to spherical aberrations, if the focus is at 
the apparent position of the molecule (we define as 2; = 0), the fiuorophore's 
defocused intensity profile's SD vs. z relation is asymmetric with respect to 
z = p. Figure S3A shows the geometry of our setup, where the direction 
of z is positive towards the glass coverslip, opposite to that of Ref. [7j. In 
our prism-type-TIRF imaging setup, the eGFP molecules were adsorbed on 
or diffusing near the fused-silica surface in TBE buffer 10.5 /xm from the 
coverslip-water interface. We used fused-silica as the TIRF interface because 
of its low background noise and thus high SNR for the study when comparing 
to objective-TIRF imaging that uses glass coverslips as the TIRF interface, 
where the mean background noise level is at least 6 times higher than ours 
at comparable laser intensity and exposure times (data not shown). 

In order to obtain accurate defocusing eGFP intensity profile parameters 
for SD and D^ calculations, we performed calculations and measurements 
of eGFP adsorbed on fused-silica surfaces. In calculations, we obtained the 
defocused fiuorophores' PSF using the diffraction integral analysis in Ref [7], 
which has been used by other groups for mainly calculating the actual axial 
location of the imaged fiuorophores [SI |9]). We assume the final PSF of a 
defocused fiuorophore to be the average of 4 emission polarizations at 0, 7r/4, 
7r/2, and 37r/4, and the light intensity at the spherical wavefront of the point 
light emitter before reaching the objective was homogeneous. Figures S3B 
and C show the SD and amplitude of the calculated fiuorophore PSFs using 
our imaging system's parameter of N.A. = 1.49, water's refractive index of 
Hi = 1.34, glass' refractive index of ^2 = 1.515, and an emission wavelength 
of 525 nm (plotted in blue). 

The experimental measurements of defocused eGFP molecules were per- 
formed with the molecules adsorbed on fused silica surfaces and a focus-drive 
(H122, Prior Scientific Inc., Rockland, MA) moving one-way in 100 nm incre- 
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ments through the focal point. The average SD and normahzed amphtude 
of eGFP intensity profiles are plotted in Figs. [7b and [Tp (in red). We 
used two fitting protocols for these results. For SD, below z = 100 nm, 
Sx,y = Sqa/I + {z/990.3y where Sq = 108.2 nm is the minimum in the eGFP 
Sx,y vs z curve that defines the focus; above z = 100 nm, a linear fit yields 
a slope of 0.73. The shape of the experimental defocusing eGFP s^^y vs z 
curve is consistent with that of the theoretical results, but the values are 
30 nm higher. These higher experimental SD values are consistent with the 
reported values in recent publications using similar imaging setups to ours 
pUl 121 [m [12], and it is due to a combination of the pixelation effect of 
the camera, finite bandwidth of the emission filter, inhomogeneity of the 
molecule's emission polarization, and the imperfection of the current single- 
molecule imaging systems. Because of this increase in the experimental eGFP 
SD values, the corresponding experimental eGFP PSF amplitudes are lower 
than the theoretical values. Below z = 150 nm, C{z) = T—T^Tiwy^] above z 

-'- + (. 726.7 ) 

= 150 nm, C{z) = . ^j^p ^ . The peak of C{z) does not coincide with the 

"'"^ 389.4 •' 

SD minimum at z = 0; rather it is shifted to -140 nm. The experimental 
eGFP functions were used for theoretical and simulation diffusing eGFP SD 
studies in this article. 

Mean emitted photon counts at each simulation step. 

The number of photons emitted at each simulation step is a random number 
drawn from a Poisson distribution with a mean value being Aexp{—z/d)C{z)s{z) 
where exp{—z/d) describes the exponentially decaying evanescent light in- 
tensity and d ^ 117 nm is the penetration depth calculated according to 
our incident angle of 70° [13], -^ is the amplitude of the refractive-index- 
mismatch affected eGFP PSF, 2ns{z)'^ corrects for the amplitude ( 2^/(^)2 )of 
a simulated Gaussian PSF with SD s{z), and A is a scaling factor that ac- 
counts for the quantum efficiency of eGFP molecules. Below we obtain C{z) 
and A. 

We first describe how C{z) is used. When we simulate a PSF by dis- 
tributing A^ photons following a 2D Gaussian spatial distribution with SD 
s{z), the amplitude of the PSF will be ^^5^)2 , where ^ ^^^ is the amplitude 
of a normalized Gaussian function with SD s{z) for one photon. However, 
the Gaussian PSFs with A^ photons in the refractive-index-mismatch case 
have the same SD s(z) but a different amplitude, ]^^ . Thus, when simu- 
lating the refractive-index-mismatch affected PSFs by spatially distributing 
photons using a 2D Gaussian distribution with SD s{z), each photon count 
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should be corrected by factor C{z)s{z)'^, where s{zy cancels the amplitude 
of the simulated normalized PSF. 

To determine A, a random value was picked to simulate PSF photon 
distributions at a finite exposure time. After the photon to camera count 
conversion (with a conversion factor M = 1 which introduces additional vari- 
ance to the emission photon distributions [2]), the modified emission photon 
count distributions were compared to the experimental distributions at the 
same exposure times. A was obtained when a good match between the two 



distributions was achieved. Figure S 10 compares the 0.6 ms experimental 



and simulated photon emission distributions; A was 0.80 and remains ap- 
proximately constant for all exposure times (the mean A value for the 0.3 
ms to 1 ms exposure times is 0.86; data not shown). 

In the theoretical Sx calculations using Eq. S|6} A was not included since 
it does not affect the final calculated PSF Sx results. 

Starting locations of the imaged diffusing eGFP molecules. 

In order to correctly simulate diffusing eGFP molecules near fused-silica sur- 
faces, the axial starting positions are needed. We obtained the eGFP diffusion 
starting position probability distributions at different exposure times by sim- 
ulating a fluorophore's emitted photon distributions for a range of starting 
positions. 

At each exposure time, we simulated 1000 axial-direction diffusion trajec- 
tories starting from the fused-silica-water interface to an extended distance 
in water for the exposure time {z = nm to 117 + 3\^2lJ^ nm measuring 
from the fused-silica surface at focus). A reflective fused-silica surface at z = 
was used. The simulations included the triplet state effect, and the num- 
ber of photons emitted at each step is described above with the difference 
of using the mean of the Poisson emission photon count distribution at each 
step, rather than drawing a random number from the Poisson distribution. 
At each starting position, we obtained the ratio of the number of photons 
emitted within the penetration depth {d = 117 nm) to all emitted photons 
for a simulated trajectory, and then the mean ratio for all simulated 1000 
trajectories was obtained and plotted in Fig. ^ The exposure times shown 
are 0.3, 0.7, and 1 ms, and half-Gaussian functions are fit to the distributions. 
The fitted SD values of the starting position distribution functions increase 
with the exposure time as SD{t) = 1.538 x lO^t + 122.4 nm. For t = 0.3 
ms. Fig. S|8]A shows that most molecules we observed experimentally should 
start within 200 nm of the surface. 
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Axial-direction PWDF^. 

In order to obtain the eGFP PWDF in the axial direction, we performed 
axial-direction diffusion simulations for all exposure times using the starting 
position distributions described above and a reflective fused-silica surface at 
z = 0. Figure S[9]A. shows 9 representative simulated PWDF^s for the 0.6 ms 
exposure time. Since 84.5% of the data can be fit by a Gaussian function 
with i?^ > 0.7, we use Gaussian function to approximate PWDF^s in Eq. q6] 
Figure S[9)3 shows the SD distribution of the fitted PWDF^s and the Gaus- 
sian fit to the distribution; the mean SD = ^/A^^QJ^si is 75.8 nm, yielding 
Az = 0.054. Az remains constant for other exposure times with a mean value 
of 0.052. Figure ^9p shows the mean value (^o) distribution of the fitted 
PWDF^s and the Gaussian fit to the distribution; the mean zo is 142.7 nm. 
The inset in Fig. SjOp shows that (zq) increases with t as {zq) = y/0. 270^1 + 
25.5 nm. For each exposure time, 1000 trajectories were simulated to obtain 
the results. 
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Figure 5: Study of the eGFP lateral PWDF^^s and their convolution with 
PSFs. (A) Nine random eGFP PWDF^^s at 0.6 ms exposure time and Gaus- 
sian fits to the unimodal distributions with R^ > 0.8. (B) The distribution 
of 1000 PWDF^ SDs, fitted with a Gaussian. (C) The 9 PWDF^s in (A) 
convolved with eGFP PSFs at focus with sq = 108.2 nm. (D) The SD dis- 
tribution of 1000 PWDF^. convolved eGFP PSFs at focus and its Gaussian 
fit. 
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Figure 6: Diffusing eGFP PCS autocorrelation plot. The black curve is a fit 
to the raw data (red dots). 
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Figure 7: (A) Our imaging setup and the schematics of emission from a 
fluorophore in water located at the fused-sihca-water interface 10.5 /zm away 
from the covershp- water interface. Dashed hues trace the emission from the 
actual location of the fluorophore, and solid lines trace the emission from the 
apparent location of the fluorophore. The letter "W" labels the wavefront 
of the emission before reaching the objective. (B) Calculated (blue) and 
experimental (red) eGFP PSF SDs^and (C) amplitude vs the defocusing 
distance z plots. Lines are flts to the experimental measurements. The focal 
point is at the minimum of the eGFP Sx vs z curve (same for Sy). 
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Figure 8: Simulation results for the diffusing eGFPs' starting-location dis- 
tribution near the fused-silica-water interface at exposure times of 0.3, 0.7, 
and 1ms (A, B, and C) and their corresponding fits. 
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Figure 9: Study of the axial-direction PWDF parameters. (A) Nine randomly 
selected PWDF^s at t = 0.6 ms. (B) Fitted-PWDF^s' SD distribution and 
its Gaussian fit. The mean is 75.8 nm. (C) Fitted-PWDF^s' mean (zq) 
distribution and the Gaussian fit. Inset, (^o) increases with Dst. 
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Figure 10: Comparing the experimental (black) and simulated (empty) pho- 
ton emission distributions at 0.6 ms. Their respective Gaussian fits in solid 
and dashed lines are in good agreement. 
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Figure 11: Comparing the experimental (green) and simulated (red) diffusing 
eGFP SD distributions. 
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